Catalytic quantum error correction 



Todd Brun* Igor Devetakt and Min-Hsiu Hsieh* 

Department of Electrical Engineering-Systems, University of Southern California, 

Los Angeles, CA 90089, USA " 

February 1, 2008 



Abstract 

We develop the theory of entanglement- asststed quantum error correcting (EAQEC) codes, 
a generalization of the stabilizer formalism to the setting in which the sender and receiver 
have access to pre-shared entanglement. Conventional stabilizer codes are equivalent to dual- 
containing symplectic codes. In contrast, EAQEC codes do not require the dual-containing 
condition, which greatly simplifies their construction. We show how any quaternary classical 
code can be made into a EAQEC code. In particular, efficient modern codes, like LDPC codes, 
which attain the Shannon capacity, can be made into EAQEC codes attaining the hashing 
bound. In a quantum computation setting, EAQEC codes give rise to catalytic quantum codes 
which maintain a region of inherited noiseless qubits. We also give an alternative construction 
of EAQEC codes by making classical entanglement assisted codes coherent. 

Information theory and the theory of error-correcting codes (coding theory) are intimately 
connected. Both address the problem of sending information over noisy channels. The sender 
Alice encodes her message as a codeword, sends it through the channel, and the receiver Bob tries 
to infer the intended message based on the channel output. 

Information theory (or rather the subfield of Shannon theory) deals with the asymptotic setting 
of increasingly long codes, with asymptotically vanishing error probability. The noisy channel 
is typically assumed to act independently on the codeword bits. The fundamental quantity of 
interest is the capacity of the channel: the optimal rate (in bits per channel use) of information 
transfer. Claude Shannon |30| gave a remarkable characterization of the channel capacity in terms 
of mutual information. Unfortunately, the capacity is achieved by random coding, which means 
highly inefficient encoding and decoding algorithms. 

Coding theory deals with the practical finite setting, characterized by a fixed code length, num- 
ber of encoded bits and correctable error set. The most popular codes have simple mathematical 
properties, such as linearity (a linear combination of codewords is another codeword), which al- 
lows for efficient encoding. The performance of these codes is then measured against the optimal 
performance set by Shannon theory. 

This relationship carries over to quantum information processing. The basic communication 
task is sending quantum information over noisy quantum channels. This setting is also relevant for 
fault tolerant quantum computation, because decoherence can be regarded as a quantum channel 
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Figure 1: The male side of the family tree of quantum Shannon theory |13| . 



connecting two points in time (rather than space). The first quantum error-correcting (QEC) code 
was discovered by Shor [STJ, leading to an explosion of research in subsequent years |TC1 [21 WR I2T1 
EsOEliniEIIB]. Calderbank and Shor [5j and Steane [33] gave the first systematic way to construct 
quantum "CSS" codes from dual-containing classical codes over Z2. These efforts culminated in 
a general theory of linear quantum codes, also known as stabilizer codes [SJEIEHIEHI- Stabilizer 
codes are equivalent to classical codes which are dual-containing with respect to the symplectic 
bilinear form. These in turn may be constructed from dual-containing classical codes over F4, 
generalizing the CSS construction [7j. 

In |31| Shor also raised the information theoretical question of characterizing the capacity of 
a quantum channel for sending quantum information, subsequently answered by |23l 1321 lll| in 
terms of coherent information. It comes as no surprise that coding theory and information theory 
continue to inform each other in the quantum setting. The capacity-achieving quantum codes of 
1 1 1 j have a structure akin to CSS codes (thanks to their common connection to cryptography). 
Concatenated stabilizer codes achieve rates equal to the coherent information evaluated on density 
operators corresponding to maximally mixed qubit states encoded by a stablizer code [21 118| . 

Research has since taken us beyond this most obvious quantum communication setting. Apart 
from quantum communication channels, there are other resources to consider, such as entanglement 
and classical communication. Great progress has been made in characterizing optimal tradeoffs 
between these resources. For example, the capacity of a quantum channel for sending classical infor- 
mation assisted by entanglement (EA capacity) is a simple single letter expression involving quan- 
tum mutual information [3J. In |13| a remarkable duality was discovered between entanglement- 
assisted quantum communication (the "father" protocol) and quantum-communication-assisted 
entanglement distillation (the "mother" protocol) . The two were shown to generate a whole family 
of protocols when combined with the more elementary protocols of superdense coding 0j , quantum 
teleportation [T] and entanglement distribution 

The father side of the family is shown in Figure ^ Quantum capacity-achieving protocols can 
be obtained from the father protocol by combining it with entanglement distribution. In conjunc- 
tion with superdense coding, the father protocol gives rise to EA capacity-achieving protocols. 
Moreover, the latter can be made coherent jl9llllllT3lll2j to recover the father protocol. 
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Can we reproduce the family in the finite setting of coding theory? Is it beneficial to do so? In 
this paper we give an affirmative answer to these two questions. We develop a general theory of 
linear "father" codes or entanglement-assisted quantum error-correcting (EAQEC) codes. The first 
and only such code to date was constructed by Bowen [H] from the [[5, 1, 3]] QEC code. EAQEC 
codes turn out to be a rather natural generalization of the usual stabilizer codes, equivalent to 
classical symplectic codes. These codes need not be dual-containing: the degree to which they are 
not measures the required amount of entanglement assistance. Consequently, any classical code 
can be made into a EAQEC code. This provides a drastic simplification, allowing the classical 
theory of error correction to be imported wholesale. 

The paper is organized as follows. Section 1 provides background on the Pauli group and 
symplectic algebra. It also reviews basic quantum strategies for sending classical information. 
Section 2 defines EAQEC codes and determines the set of errors they can correct. Section 3 
generalizes the code construction method of |Hj based on classical codes over F4. Section 4 
regards the right branch of Figure ^ constructing catalytic QEC codes from EAQEC codes. 
Section 5 regards the left branch of Figure ^ constructing entanglement assisted codes for sending 
classical information (EACEC codes). These are then made coherent providing an alternative 
construction of EAQEC codes. Section 6 discusses bounds on the performance of EAQEC codes. 
Section 7 recovers Bowen's result in our framework. Section 8 updates the table of known codes 
from jS). We discuss our results in Section 9. 

1 Background 

In this section we review the properties of Pauli matrices, and relate them to symplectic binary 
and quaternary vector spaces. Our presentation follows Forney et al. [23] and Hamada |I8| . 

1.1 Single qubit Pauli group 

A qubit is a quantum system corresponding to a two dimensional complex Hilbert space Ti.. Fixing 
a basis for H,, the set n of Pauli matrices is defined as 
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The Pauli matrices are Hermitian unitary matrices with eigenvalues belonging to the set {1,-1}. 
The multiplication table of these matrices is given by: 
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Observe that the Pauli matrices either commute or anticommute. Let [A] = {/3A \ f3 € C, \ j3\ = 
1} be the equivalence class of matrices equal to A up to a phase factor. 1 Then the set [II] = 
{[I], [X], [Y], [Z]} is readily seen to form a commutative group under the multiplication operation 
defined by [A][.B] = [AB]. It is called the Pauli group. 

We are interested in relating the Pauli group to the additive group (Z2) 2 = {00, 01, 10, 11} of 

1 It makes good physical sense to neglect this overall phase, which has no observable consequence. 
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binary words of length 2 described by the table: 
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This group is also a two-dimensional vector space over the field %i- A bilinear form can be defined 
over this vector space, called the symplectic form or symplectic product 2 : (Z2) 2 x (Z2) 2 — > 1i, 
given by the table 
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In what follows we will often write elements of (Z2) 2 as u — {z\x), with z, x € Z2. For instance, 01 
becomes (0|1). For u = (z\x),v = {z'\x') G (Z2) 2 the symplectic product is equivalently defined by 

u © v — zx' — z'x. 

Define the map N : (Z 2 ) 2 — > II by the following table: 
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This map is defined in such a way that N( z \ x ) and Z Z X X are equal up to a phase factor, i.e. 

[N {z]x) ] = [Z*X% 

We make two key observations 

1. The map [N] : (Z 2 ) 2 — > [II] induced by N is an isomorphism: 

[N U ][N V ] = [N u+V \. 

2. The commutation relations of the Pauli matrices are captured by the symplectic product 

N U N V = {~l) uQv N v N u . 

Both properties are readily verified from the tables. 
1.2 Multi-qubit Pauli group 

Consider an n-qubit system corresponding to the tensor product Hilbert space TL® n . Define an 
n-qubit Pauli matrix A to be of the form A = A\ ® A2 ® • • • A n , where Aj e II. The set of all 
4™ n-qubit Pauli matrices is denoted by II". The product of elements of II™ is an element of II™ 
up to a phase factor. Define as before the equivalence class [A] = {/3A | (3 £ C, = 1}. Then 

[A][B] = [Aiffi ] © [A 2 B 2 ] ® ■ ■ ■ ® [A n B n ] - [AB]. 

2 Strictly speaking it is not an inner product. 
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Thus the set [II"] = {[A] : A £ 11"} is a commutative multiplicative group. 

Now consider the group/vector space (Z2) 2 " of binary vectors of length 2n. Its elements may 
be written as u = (z|x), z = Z\ . . . z n £ (Z2)™, x = x\ . . . x n £ (Z2)". We shall think of u, z and x 
as row vectors. The symplectic product of u — (z|x) and v = (z'|x') is given by 

„ T iT I T 

u V = zx - Z X . 

The right hand side are binary inner products and T denotes the transpose. This should be thought 
of as a kind of matrix multiplication of a row vector and a column vector. We use u v T rather 
than the more standard u v T to emphasize that the symplectic form is used rather than the binary 
inner product. Equivalently, 

U V T = ^ Uj V j 

i 

where ttj = (zi\xi), Vi = (z'^x^) and this sum represents Boolean addition. Observe that u0u T = 0, 
i.e., every vector is "orthogonal" to itself. 

The map N : (Z2) 2 ™ — » IT" is now defined as 

N u = N Ul ®-~®N Un . 

Writing 

X x = X Xl ®---®X x ", 
Z z = Z Z1 

as in the single qubit case, we have 

[iV (z | x) ] = [Z Z X*]. 
The two observations made for the single qubit case also hold: 

1. The map [N] : (Z2) 2 " — > [II"] induced by N is an isomorphism: 

[iV u ][iV v ] = [iV u+v ]. (1) 

Consequently, if {ui, . . . ,u m } is a linearly independent set then the elements of the Pauli 
group subset {[iV Ul ], . . . , [-/V Um ]} are independent in the sense that no element can be written 
as a product of others. 

2. The commutation relations of the n-qubit Pauli matrices are captured by the symplectic 
product 

N U N V = (-l) uQvT N V N U . (2) 



1.3 Properties of the symplectic form 

In this subsection we present two results which will play a major role in the construction of 
EAQEC codes. Together they will enable us to conclude that any independent subset of the n- 
qubit Pauli group can be transformed via a unitary operation into a canonical set whose elements 
act nontrivially only on single qubits. 

A subspace V of (Z2) 2 " is called symplectic ^Uj if there is no v £ V such that 

v © u T = 0, Vu £ V. (3) 
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(Z2) 2 ™ is itself a symplectic subspace. Consider the standard basis for (Z2) 2 ™, consisting of gi = 
(ei|0) and h; = (0|ej) for i = 1, . . . , n, where = (0, . . . , 0, 1, 0, . . . , 0) [1 in the ith position] are 
the standard basis vectors of (Z2)™. Observe that 

gi©gj = 0, forain,j (4) 

hi©hj = 0, foralli,j (5) 

gi 0hJ = O, foraUi^j (6) 

g, hf = 1, for all i. (7) 

Thus, the basis vectors come in n hyperbolic pairs (gj,hj) such that only the symplectic product 
between hyperbolic partners is nonzero. The matrix J — [gj hj] defining the symplectic product 
with respect to this basis is given by 

j / Onxn InXn 

\ Inxn 0„ x n 

where I nxn and nxn are the n x n identity and zero matrices, respectively. A basis for (Z2) 2 " 
whose symplectic product matrix J is given by (JSJ is called a symplectic basis. In the Pauli picture, 
the hyperbolic pairs (gj, hj) correspond to (Z ei , X ei ) - the anticommuting Z and X Pauli matrices 
acting on the ith qubit. 

In contrast, a subspace V of (Z2) 2n is called isotropic if J2J holds for all v e V. The largest 
isotropic subspace of (Z2) 2 ™ is n-dimensional. The span of the gi, i = 1, . . . , n, is an example of a 
subspace saturating this bound. 

A general subspace of (Z2) 2 ™ is neither symplectic nor isotropic. The following theorem, stated 
in PI]] and rediscovered in Pauli language in \T^, says that an arbitrary subspace V can be 
decomposed as a direct sum of a symplectic part and an isotropic part. We give an independent 
proof here. 

Theorem 1.1 Let V be an m- dimensional subspace of (Z2) 2 ™ . Then there exists a symplectic basis 
of (Z2) 2 ™ consisting of hyperbolic pairs (uj, Vj), i = 1, . . . , n, such that {ui, . . . , u c+ £, Vi, . . . , v c } 
is a basis for V , for some c, t > with 2c + £ = m. 
Equivalently, 

V = symp(y) iso(V) 

where symp(V r ) = spanjui, . . . , u c , vi, . . . , v c } is symplectic and iso(V) = span{u c+ i, . . . , u c+ i} is 
isotropic. 

Remark It is readily seen that the space iso(K) is unique, given V. In contrast, symp(V^) is not. 
For instance, replacing vi by v' x = vi + u c+ i in the above definition of symp(V) does not change 
its symplectic property. 

Proof Pick an arbitrary basis {wi, . . . , w TO } for V and extend it to a basis {wi, . . . , W2„} for 
(Z2) 2n . We will describe an algorithm which yields the basis from the statement of the theorem. 

The procedure consists of n rounds. In each round a new hyperbolic pair (uj,Vj) is generated; 
the index i is added to the set U (V) if Ui e V (v.; 6 V) . 

Initially set i = 1, m' = m, and U = V = 0. The ith round reads as follows. 

1. We start with vectors wi, . . . , w 2 ( n _, + i), and m, . . . Uj-x, vi, . . . Vj_i, such that 

(a) wi, . . . , w 2 („_. i+ i), ui, . . .Uj_i,vi, . . .Vj_i is a basis for (Z 2 ) 2 ™, 

(b) each of Ui, . . . Uj_i,Vi, . . . v»_i has vanishing symplectic product with each of Wi, . . . , w 2 („_ i+ i), 

(c) V = spanjwj : 1 < j < m'} spanjuj :j6W}© spanjvj : j e V}. 
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These conditions are satisfied for i = 1. 

2. Define u.; =w%. If to' > I then and add i to U. Let j > 2 be the smallest index for which 
Wi wj = 1. Such a j exists because of (a), (b) and the fact that there exists a w <s (Z 2 ) 2n 
such that Ui w T = 1. 

Set v, = wj. 

3. If j < to': 

This means that there is a hyperbolic partner of in V. Add i to V; swap Wj with W2; for 
k = 3, . . . , 2{n — i + 1) perform 

w fc-2 : = W A; - (Vi W^)U, - (Ui © Wfc )Vi, 

so that 

w' fe _2 © uf = w' fe _ 2 vf = 0; (9) 

set to' := to' — 2. 
If j > to': 

This means that there is no hyperbolic partner of in V. Swap Wj with w 2 („_ i+1 ); for 
k = 2, . . . , 2(n — i) + 1 perform 

W 'fc-1 ; = W fc - ( V * © W fc ) U » - ( U * © W fc ) V »: 

so that 

w' fe _ 10 uf =w' fe _ 10 vf = 0; (10) 
if to' > 1 then set m' := to' — 1. 

4. Let Wfc := wl for 1 < fc < 2(n — i). We need to show that the conditions from item 1 are 
satisfied for the next round (i := i + 1). Condition (a) holds because {iij, Vi,w[, . . . wL n _ i \} 
are related to the old {wi, . . . W2( n _j+i)} by an invertible linear transformation. Condition 
(b) follows from © and l |10|l . Regarding condition (c), if to' = then it holds because U 
and V did not change from the previous round. Otherwise, consider the two cases in item 3. 
If j < to' then {iij, Vj, w' 1; . . . w' Tl ,_ 2 } are related to the old {wi, . . . w m /} by an invertible 
linear transformation. If j > to' then {iij, w' 1; . . . w^ n ,_ 1 } are related to the old {wi, . . . w m /} 
by an invertible linear transformation (the (iij © w^)vi terms vanish for 1 < k < m! because 
there is no hyperbolic partner of in V). 

< to' < 2(n — i) at the end of the ith round. Thus to' = after n rounds and hence V = 
spanjuj : j £ U} © spanjvj : j £ V}. The theorem follows by suitably reordering the (iij-, Vj). □ 

A symplectomorphism T : (Z2) 2 ™ — > (Z2) 2n is a linear isomorphism which preserves the sym- 
plectic form, namely 

T(u)©T(v) T = u©v T . (11) 

The following theorem relates symplectomorphisms on (Z 2 ) 2 " to unitary maps on 7i® n . It 
appears, for instance, in For completeness, we give an independent proof here. 

Theorem 1.2 For any symplectomorphism T on (Z2) 2n there exists a unitary map Ur on 7i®" 
such that for all u £ (Z 2 ) 2 ™, 

[iV T(u) ] = [C/xA^u^ 1 ]. 
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Remark. The unitary map Ur may be viewed as a map on [IT] given by [A] i— > \UyAU T 1 ]. The 
theorem says that the following diagram commutes 



(Z 2 ) 2 " 

[N] 

[n] 



(Z 2 ) 2 « 
[ [N] 

[n] 



Proof Consider the standard basis gi = (ej|0), = (0|ej). Define the unique (up to a phase 
factor) state |0) on 7i®" to be the simultaneous +1 cigcnstate of the commuting operators N gj , 
j = 1, . . . ,n. Define an orthonormal basis {|b) : b = b\ . . . b n G (Z2)™} for 7i®" by 

|b>=AT EjMli |0>. 

The orthonormality follows from the observation that |b) is a simultaneous eigenstate of N gj , 
j = 1, ... ,71 with respective eigenvalues (— l) bi : 

N Sj \b)=N gj N Eihhi \0) 

= (-lpN^ tb ^N S] \0) 
= (-l)^iV EjMlj |0) 
= (-l) 6 '|b). 



(12) 



The second line is an application of (J2J. 

Define gj := T(gj). We repeat the above construction for this new basis. Define the unique 
(up to a phase factor) state |0) to be the simultaneous +1 eigenstate of the commuting operators 
N Si , i — 1, . . . , 71. Define an orthonormal basis {|b)} by 



l b > =Ar £ j6i hJ >- 

Defining u = ^ z iSi + x i^i, u = J2i z i&i + x i^i and x = x i ■ ■ ■ x n, we have 



(13) 



(14) 



AT a |b)=iV a iV Eibii;i |0) 

= ( _ 1) a0(E i a 4 r A r E ^ ii . iiV ~|o) 
= (_i)ao(E I fc.h 1 ) T e ie(S) 7Vj 

= (_i)S®<E«M^) T e <fl(s)jv_ 

= (_i)S©(E«fc£*) T e *(5)|b^) 

= (_i)»o(E«M^ T e *(s)|b^), 

where 6*(u) is a phase factor which is independent of b. The first equality follows from 

the second from (J2J, the third from QJ, the fourth from the definition of |0) and the fact that 

X b X x = X h+X , the fifth from O, and the sixth from (JTTJ- Similarly 



N u \b) = (-l) u0 ^. b ' h ') T e iv(u) |b + x), 

where v?( u ) is a is a phase factor which is independent of b. 
Define Ur by the change of basis 

(y T = ]T|b)<b|. 



(15) 
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Combining 114(1 and l|15[) gives for all |b) 

N r{u) U r \b) = (-l) u0 E. b ^) T e ie(il) C/T|b + x) 
= e i[e(U) -^ u)] U r N u \b). 

Therefore [JV T (u)] = [UyNnUy 1 ]. 

1.4 Encoding classical information into quantum states 

In this subsection we review two schemes for sending classical information over quantum channels: 
elementary coding and superdense coding. These will be used later in the context of quantum 
error correction to convey information to the decoder about which error happened. 

In the first scheme, Alice and Bob are connected by a perfect qubit channel. Alice can send an 
arbitrary bit a £ Z2 over the qubit channel in the following way: 

• Alice locally prepares a state |0) in TC. This state is the +1 eigenstate of the Z operator. 
Based on her message a, she performs the encoding operation X a , producing the state \a). 

• Alice sends the encoded state to Bob through the qubit channel. 

• Bob decodes by performing the von Neumann measurement in the {|0), |1)} basis. As this 
is the unique eigenbasis of the Z operator, this is equivalently called "measuring the Z 
observable" . 

We call this protocol "elementary coding" and write it symbolically as a resource inequality |12l 

nam 3 

[«->• q] > [c-> c]. 

Here [q — ► q] represents a perfect qubit channel and [c — » c] represents a perfect classical bit 
channel. The inequality > signifies that the resource on the left hand side can be used in a 
protocol to simulate the resource on the right hand side. 

Elementary coding immediately extends to m qubits. Alice prepares the simultaneous +1 
eigenstate of the Z ei , . . . , Z Bm operators |0), and encodes the message a S {1>2) m by applying A a . 
Bob decodes by simultaneously measuring the Z ei , . . . , Z Gm observables. We could symbolically 
represent this protocol by 

m [q — > q] > m [c — ► c] . 
In the second scheme, Alice and Bob share the ebit state 

|#) = -^(|0)®|0) + |1)«|1)) (17) 

in addition to being connected by the qubit channel. In l|17l) Alice's state is to the left and Bob's 
is to the right of the ® symbol. 

The state |$) is the simultaneous (+1,+1) eigenstate of the commuting operators Z ® Z and 
X®X. Again, the operator to the left of the (£> symbol acts on Alice's system and the operator to 
the right of the ® symbol acts on Bob's system. Alice can send a two-bit message (ax, aq) S (Z2) 2 
to Bob using "superdense coding" 0]: 

• Based on her message (01, 02), Alice performs the encoding operation Z ai X a2 on her part of 
the state |$). 

3 In II 21 resource inequalities were used in the asymptotic sense. Here they refer to finite protocols, and are thus 
slightly abusing their original intent. 
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Figure 2: A generic entanglement assisted quantum code. 

• Alice sends her part of the encoded state to Bob through the perfect qubit channel. 

• Bob decodes by performing the von Neumann measurement in the {(Z ai X a2 (g J)|$) : 
(01,^2) <= (2<2) 2 } basis, i.e., by simultaneously measuring the Z ® Z and A <g A observ- 
ablcs. 



The protocol is represented by the resource inequality 

[q^q] + [qq]>2[c 



(18) 



where [qq] now represents the shared ebit. It can also be extended to m copies. Alice and Bob 
share the state |<!>)® m which is the simultaneous +1 eigenstate of the Z ei (g Z ei , . . . , Z Bm (g Z Bm 
and X ei <g X ei , . . . , A 6m (g A 6m operators. Alice encodes the message (ai,a 2 ) € (Z 2 ) 2m by 
applying A ai Z &2 . Bob decodes by simultaneously measuring the Z ei (g) Z ei , . . . , Z Bm (g) Z 6m and 
A ei (g A ei , . . . , X Bm (g A 6m observables. The corresponding resource inequality is 

m[q —> q] + m[qq] > 2m [c — > c] . 

Superdense coding provides the simplest illustration of how entanglement can increase the 
power of information processing. 



2 Entanglement-assisted quantum error correction 

In this section we introduce entanglement-assisted error correcting (EAEC) codes and prove our 
main result, Theorem 2.4, which gives sufficient error-correcting conditions. 

2.1 The model: discretization of errors 

It is well known that for standard quantum error correction (i.e., that unassisted by entanglement) 
it suffices to consider errors from the Pauli group (see e.g. [2HI-) We will show this for entanglement- 
assisted quantum error correction. Denote by C the space of linear operators defined on the 
qubit Hilbert space Tt. We will often encounter isometric operators U : — > The 
corresponding superoperator, or completely positive, trace preserving (CPTP) map, is marked by 
a hat U : C® ni -> £®" 2 and defined by 

U(p) = UpU ] . 
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Observe that U is independent of any phases factors multiplying U. Thus, for a Pauli operator 
N u , N u only depends on the equivalence class [N u ]. 

Our communication scenario involves two spatially separated parties, Alice and Bob, as depicted 
in Figure [3 The resources at their disposal are 

• a noisy channel defined by a CPTP map Af : C® n — > C® n taking density operators on Alice's 
system to density operators on Bob's system; 

• the c ebit state |$)® c shared between Alice and Bob. 

Alice wishes to send k qubits perfectly to Bob using the above resources. An [[n, k; c]] entanglement- 
assisted quantum error correcting (EAQEC) code consists of 

• An encoding isometry E = U cnc : C® k ® C® c — > C® n 

• A decoding CPTP map V : £®« <g> C® c -> £® k 
such that 

VoM o U cnc o U app = id® fe , 
where U app is the isometry which appends the state |<!>)® c , 

c/ apP M = |^)|$)® c , 

and id : C — ► C is the identity map on a single qubit. The protocol thus uses up c ebits of 
entanglement and generates k perfect qubit channels. We represent it by the resource inequality 
(with a slight abuse of notation \V2\) 

(Af) +c[qq] >k[q^q}. 

Even though a qubit channel is a strictly stronger resource than its static analogue, an ebit of 
entanglement, the parameter k — c is still a good (albeit pessimistic) measure of the net noiseless 
quantum resources gained. It should be borne in mind that a negative value of k — c still refers to 
a non-trivial protocol. 

To make contact with classical error correction it is necessary to discretize the errors. This is 
done in two steps. First, the CPTP map Af may be (non-uniquely) written in terms of its Kraus 
representation 

i 

Second, each Aj may be expanded in the Pauli operators 

ue(z 2 ) 2 " 

Define the support of Af by supp(AT) = {u € ^Li) 2n : 3i, ai yU ^ 0}. The following theorem allows 
us, absorbing U app into U cnc , to replace the continuous map Af by the error set S — supp(A/"). 

Theorem 2.1 IfVoN u o U cnc = id® k for all u G supp(A/") 7 then V o Af a U cnc = id® fe . 

Proof We may extend the map T> to its Stinespring dilation - an isometric map C/dec with a 
larger target Hilbert space C® k ® such that 

V = Tic ot) dcc . 

The premise of the theorem is equivalent to saying that for all u £ supp(A/") and all pure states 
\<p) in H® n , 

Ud cc N u U cnc \tp) = \<p) O |u) 
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for some pure state |u)(u| on By linearity 

Ud ec AiU enc \(p) = \<p) <g> \i), 
with the unnormalized state |i) = X)u a *,ul u )- Furthermore, 

{Udoc oTVo U enc )(\ip)(tp\) = U dcc (^AiU^MVLM) U L 

\ i I (19) 

= \<pM ® 

where the second subsystem corresponds to Tracing out the latter gives 

(VoMoU cnc )(\^\) = |^|, 
concluding the proof. □ 



2.2 The canonical error set and syndrome coding 

By the results of the previous subsection, we are now interested in EAQEC codes which correct a 
particular error set S C (Z 2 ) 2 ". We first restrict attention to a simple error set, which will turn 
out to be generic due to the results of Section IT31 

Consider the following trivial encoding operation Uq defined by 

C/o:|^>|$)® c ^|^)|0)|a>)^. (20) 

In other words, the register containing |0) (of size £ = n — k — c qubits) is appended to the registers 
containing \<p) (of size k qubits) and |<!>)® c (of size c qubits each for Alice and Bob). What errors 
can she correct with such a simple-minded encoding? 



Proposition 2.2 The code given by Uo and a suitably defined decoding map T>q can correct the 
error set 

S = {(a(a,ai,a 2 ),b,a 2 |/3(a,ai,a 2 ),a,ai) : b,a e (Z 2 )*,ai,a 2 G (Z 2 ) c }, (21) 
for any functions a, (3 : (Z 2 )^ x (Z 2 ) c x (Z 2 ) c — > (Z 2 ) fc . 

Proof The protocol is shown in figure |3| After applying an error iV u with 

u = (a(a,ai,a 2 ),b,a 2 |/?(a,ai,a 2 ),a,ai), (22) 
the encoded state |<p)|0)|<I>) l2lc becomes (up to a phase factor) 

Z a(a,a 1 ,a 2 ) x ^(a,a 1 ,a 2 )|^ ^ X a Z h \0) <g> (Z ai A &2 (g> /)|$)^ c 
= Z a(a,a ll a a )j f /J(a,a 1 ,a a )|^ g | a } g | ai ,a 2 ), 

where |a) = X a \0) and |a x ,a 2 ) = (Z ai X a2 <£> 7)|$)® c . As the vector (a,a!,a 2 ,b) T completely 
specifies the error u, it is called the error syndrome. The state (|23[) only depends on the reduced 
syndrome r = (a, ai,a 2 ) T . In effect, a and (ai,a 2 ) have been encoded using plain and superdense 
coding, respectively. Bob, who holds the entire state l|23|l . may identify the reduced syndrome using 
the results of section 11.41 Bob simultaneous measures the Z ei , . . . , Z Bm observables to decode a, 
the Z ei (Xi Z ei , . . . , Z e - <x> Z e ™ observables to decode ai , and the X ei <x> X ei , . . . , X e ™ <x> X e ™ 
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\<t>) 



\<p) 

10) 



k qubits 



t qubits 



c qubits 



\ 



x a z b 



x- a z~P- 



» 



{\a)(a\}- 



{|a 1 a 2 )<a 1 a 2 |} 



c qubits 



a\d2 



Figure 3: The canonical code. 

observables to decode a 2 . He then performs Z~ Q ( a > a i> a 2) X _ ' 3 ' a ' ai,a2 J on the remaining fc qubit 
system, leaving it in the state 

Since the goal is the transmission of quantum information, no actual measurement is necessary. 
Instead, Bob can perform the CPTP map X>o consisting of the controlled unitary 

£/odcc= J2 ^" a(a,ai,a2) ^"^ (a,ai ' a2) ®|aXa|®|a 1) a 2 )(a 1 ,a 2 |, 



followed by discarding the last two subsystems. 



□ 



The above code is degenerate with respect to the error set S, which means that the error can 
be corrected without knowing the full error syndrome. 

We can characterize our code in terms of the parity check matrix F given by 



F = 



Fi 



Fi — ( Otxk iixi Oexc | Oexk ®ixi Otxc ) , 
F s = 



Ocxk Ocxi Icxc 
O c xfc ®cx£ CXC 



Ocxk Ocxl CX c 
O c xfe CX £ Icxc 



(24) 
(25) 
(26) 



with I = n — k — c. 

The vector space rowspace(-F) decomposes into a direct sum of the isotropic subspace rowspace(-Fj) 
and symplectic subspace rowspace(_F,g), as in Theorcm ll.il Define the symplectic code correspond- 
ing to F by 

Cq = rowspace(F) 



where 

Note that (V J 
rowspace(Fs). 



= {w : w u T = 0, Vu e V}. 
V. Then Cq = rowspace(F), iso(Cg ) = rowspace(Fj) and symp(C 1 
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The number of ebits used in the code is 



c = — dimrowspace(i 7 5) 

and the number of encoded qubits is 

k = n — dimrowspace^/) — — dimrowspace^s). 

The parameter k — c which is the number of encoded qubits minus the number of ebits used is 
independent of the symplectic structure of F: 

k — c — n — dimrowspace(-F). 

The error set So can be described in terms of F: 

Proposition 2.3 The set So of errors correctable by the code So is such that, if u, u' S So and 
u/u', then either 

i) u — u' ^ Co (equivalently: F (u — u') T ^ T ), or 

ii) u — u' e \so{Cq) (equivalently: u — u' € rowspace(i 7 /) / ) 

Proof If u is given by ill' 21) then F u T = r = (a, ai,a2) T , the reduced error syndrome. 
By definition (|21|l . two distinct elements of So either have different reduced syndromes (a, ai,a2) 
(condition 1) or they differ by a vector of the form (0, b, 0|0, 0, 0) (condition 2). Observe that 
condition 1 is analogous to the usual error correcting condition for classical codes [27] • □ 

The parity check matrix F also specifies the encoding and decoding operations. The space 7^® fe 
is encoded into the codespace defined by 

C = {Uo\vm® c :\v)£n® k }. 

It is not hard to see that the codespace is the simultaneous +1 eigenspace of the commuting 
operators: 

1. I®Z e >®I(E)I, i = l,...J; 

2. I®J®Z e >8Z e >, j = l,...,c; 

3. /«/®r'8r', j = 1, . . . , c. 

Above, the first three operators act on Alice's qubits and the fourth on Bob's. Define the matrix 



B = 

Define the augmented parity check matrix 
^aug = (F, B) = 



Of xc 


Oixc 


lex c 


Ocxc 


Ocxc 


Icxc 







Ofxc 


Ofxc 


Oixk 


Oe xi 


Oixc 


Of xc 


c xfc 


Ocxf 


Icxc 


Icxc 


Ocxfc 


Ocxe 


O c xc 


O c xc 


Ocxfc 


c x£ 


O c xc 


O c xc 


Ocxfc 


O c xl 


Icxc 


Icxc 



(27) 



(28) 



Observe that rowspace(-F aug ) is purely isotropic. The codespace is now described as the simulta- 
neous +1 eigenspace of 

{A w : w S rowspace(F aug )}, 

or, equivalently that of 

Go = {A^ w : w is a row of -Fkug}- 
The decoding operation T>o is also described in terms of F. The reduced syndrome r = F0 u T is 
obtained by simultaneously measuring the observables in Go- The reduced error syndrome corre- 
sponds to a number of possible errors u G So which all have an identical effect on the codespace. 
Bob performs A u = A_ u to undo the error. 
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2.3 The general case 

We now present our main result: how to convert an arbitrary (n + fc)-dimensional subspace C of 
(Z2) 2n into a EAQEC code. Consider the (n — fc)-dimensional subspace C . By Theorem 11.11 
there exists a symplectic basis of (Z2) 2 " consisting of hyperbolic pairs (u^, v^), i = 1, ...,n, such 
that the ordered set 1Z = {uk+i, . . . , u„, Vk+e+i, ■ ■ ■ , v„} is a basis for C , for some c, I > with 
2c + £ = n — k, and k — c = k. Let H be the matrix whose rows consist of the elements of 1Z in the 
order given from top to bottom. Let T be the symplectomorphism defined by 

T(iii) = g, (29) 
T(vi) = hi. (30) 

Recall the matrix F given by (|24(I - H26|I . Observe that, with a slight abuse of notation, 

T(-ff) = F 

in the sense that T takes the ith row of H to the zth row of F. We may extend T to act on 
(Z 2 ) 2 ( n+C \ including a trivial action on the bits corresponding to Bob's side. Then 

Y(i?aug) = F^ug, (31) 

where i? aug = (H,B). 

In terms of vector spaces 

T(C X ) = C (32) 
T(iso(C ± )) = iso(C L ). (33) 

Note that c = \ dimsymp(C x ). We our now ready for our main result: 

Theorem 2.4 There exists an [[n, k;c]] EAQEC code (U enc ,T>) with the following properties: 

1. It can correct the error set S defined by: if u, u' 6 S and u/u', then either 

i) u — u' ^ C (equivalently: H (u — u') T 7^ T ), or 
ii) u — u' £ iso(C _L ) (equivalently: u u' G rowspace(ffj) J. 

2. The codespace C — C/enc(W® ) is a simultaneous eigenspace of the ordered set 

Q = {N w : w is a row of iJ aU g}, 
where H aug — (H,B), with B given by \27j) . 

3. To decode, the reduced error syndrome 

r = H u T (34) 

is obtained by simultaneously measuring the observables from Q. Bob finds a u satisfying 
\34{ ) and performs N u to undo the error. 

Remark The above theorem generalizes the error correcting conditions of ^| |H] for quantum 
error correcting codes unassisted by entanglement. When c — then C = iso(C ± ) and no 
entanglement is used in the protocol. We call such codes dual- containing. 

Proof By Theorem II . 21 there exists a unitary Uy such that for all u £ (Z2) 2 ™ 

[7V T(U) ] = [UrN^ 1 ], (35) 

and hence 

The above also holds for T and Uy extended to act trivially on Bob's side. 

Our EAQEC code is defined by U crLC = U^ 1 Uo and T> = T>o o Uy : as shown in Figure 0| 
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|0>C 















x a z fc 


■u- 1 


V 










V 










V 





x -a z -@ 



{|a 1 a 2 >(a 1 a 2 |} 



Figure 4: Generalizing the canonical code construction. 

1. Recall the error set So defined in Proposition ^. 31 From (|32[l and i|33|) it follows that T(S) = 
Sq. By Proposition 12. 31 for all u G S, 

V o A T(u) o U = id®\ 

from which 

T> o A u o L> cnc = id 8fc 
follows. Thus, the code (U enc ,T>) corrects the error set S. 

2. The codespace is C = U^ 1 (C ), by definition. According to C is the simultaneous +1 
eigenspace of 

Go = {-ZVt(w) : w is a row of -ff aug }, 

or by the set 

= {UrNv/Uy 1 : w is a row of i? a ug}- 
Lemma 12.51 now implies that the codespace C is a simultaneous eigenspace of Q . 

3. Assume that error u£5 occurs. The operation T>q involves 

(a) measuring the set of operators given by Go, or equivalently Go, yielding the reduced 
syndrome 

r = F0T(uf; 

(b) performing Ay( u ), where T(u) S Sq is an error consistent with the observed syndrome 



(I3H) holds because 



r = T(i7)OT(u) T = #0u 2 



By lemma T2.6I performing T> = Dq o Ut is equivalent to measuring the set G, followed by 
performing N u = Uy 1 o iV-r(u) ^x, followed by Ur to undo the encoding. If the final Ur is 



omitted, one recovers the encoded state rather than the original one. 



□ 



Lemma 2.5 If Co is a simultaneous eigenspace of Pauli operators from the setG'o thenC = U ^Co) 
is a simultaneous eigenspace of Pauli operators from the set G — {t/ -1 Ai7 : A € G'q}- 



lG 



Proof Observe that if 

A\il>) = a\il>), 

then 

(U^AUp-^ip) = aU~ l \i>). 

□ 

Lemma 2.6 Performing U followed by measuring the operator A is equivalent to measuring the 
operator U~ 1 AU followed by performing U . 

Proof Let IL be a projector onto the eigenspace corresponding to eigenvalue of A. Performing 
U followed by measuring the operator A is equivalent to the instrument (generalized measurement) 
given by the set of operators {ILL/}. The operator U~ 1 AU has the same eigenvalues as A, and the 
projector onto the eigenspace corresponding to eigenvalue is L/ _1 ILL/. Measuring the operator 
U~ 1 AU followed by performing U is equivalent to the instrument {L/(L7 _:L ILL/)} = {ILL/}. □ 

2.4 Distance 

The notion of distance provides a convenient way to characterize the error correcting properties 
of a code. We start by defining the weight of a vector u = (z|x) € (Z2) 211 by wt(u) = wt(z V x). 
Here V denotes the bitwise logical "or", and wt(y) is the number of non-zero bits in y G (Z2)™. In 
terms of the Pauli group, wt(u) is the number of single qubit Pauli matrices in N u not equal to 
the identity I. 

Consider a symplectic code C. The distance of C is the maximum d such that for each nonzero 
u of weight < d cither 

i) u £ C, or 

ii) u G iso(C x ) 

It is called non- degenerate if the second condition is not invoked. A code is said to correct t 
errors if it corrects the error set {u : wt(u) < t} but not {u : wt(u) < t + 1}. Comparing these 
definitions with Theorem 12.41 a code with distance d = 2t + 1 can correct t errors. An [[n, k; c]] 
EAQEC code with distance d will be referred to as an [[n, k, d; c]] code. 

3 Relation to quaternary codes 

We shall now show how to construct non-degenerate EAQEC codes from classical codes over F4, 
generalizing the work of [H]- Following the presentation of Forney et al. [23, the addition table of 
the additive group of the quaternary field F4 = {0, 1,lu,lJ} is given by 



+ 





CO 


1 


UJ 








UJ 


1 


UJ 


UJ 


~ 





UJ 


1 


1 


1 


CO 





UJ 


UJ 


CO 


1 


Zj 






Comparing the above to the addition table of (Z2) 2 establishes the isomorphism 7 : F4 — * (Z2) 2 , 
given by the table 



F 4 







00 


UJ 


01 


1 


11 


UJ 


10 



17 



The multiplication table for F4 is defined as 



X 





UJ 


1 


UJ 

















UJ 





UJ 


57 


1 


1 





57 


1 


UJ 


UJ 





1 


UJ 


UJ 



Define the traces (Tr) of the elements {0,1, uj,uj} of F4 as {0,0,1,1}, and their conjugates ("'") 
as {0, 1,uj,uj}. Intuitively, Tra measures the "w-ness" of a G F4. Observe that a = if and only if 
both Trcua = and Tr uja — 0. The Hermitian inner product of two elements a, b G F4 is defined 
as (a, b) = a^b G F4. The trace product is defined as Tr(a, b) G F2. The trace product table is 
readily found to be 



Tr(,> 





UJ 


1 


UJ 

















UJ 








1 


1 


1 





1 





1 


UJ 





1 


1 






Comparing the above to the table of (Z2) 2 establishes the identity 

Tr(a,6)= 7 (a)0 7 (6). 

These notions can be generalized to n-dimensional vector spaces over F4. Thus, for a, b G (F 4 )™, 

Tr(a,b> = 7 (a)0 7 (b) T . (36) 
Let wt4(a) be the number of non-zero bits in a G (F4)™. Then we have another identity 

wt( 7 (a)) = wt4(a), (37) 

where 7 (a) G (Z 2 ) 2 ". 

Proposition 3.1 If a classical [n,k,d]4 code exists then an [[n, 2k~ n + c, d; c]] EAQEC code exists 
for some non-negative integer c. 

Proof Consider a classical [n, k, d]i code (the subscript 4 emphasizes that the code is over F4) 
with an (n — k) x n quaternary parity check matrix H4. By definition, for each nonzero a G (F4) n 
such that wt4(a) < d, 

(H i ,a)^0 T . 

This is equivalent to the logical statement 

Tr(w# 4 ,a) ^0 T V Tr(aJJJ 4 , a) ^ T . 

This is further equivalent to 

Tr(tf 4 ,a) ^0 T , 

where 

Define the (2n — 2k) x 2n symplectic matrix H — -7(^/4). By the correspondences (|36|l and (|37|l . 

HQu T ^0 T , 
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holds for each nonzero u 6 (Z^) 2 ™ with wt(u) < d. Thus C = rowspace(i/) J - defines a non- 
degenerate [n, 2k — n + c, d; c] EAQEC code, where 

c = — dimsymp(C). 

□ 

Any classical binary [n, k, d]a code may be viewed as an quaternary [n, k, d]±. In this case, the 
above construction gives rise to a CSS-type code. 



4 Catalytic quantum error correcting codes 

So far we have been considering communication scenarios involving two spatially separated parties 
Alice and Bob connected by a noisy channel Af . In this setting, entanglement between them is a 
meaningful resource. What if Alice and Bob are separated only in time? For example, Af could 
represent the time evolution of the state of a quantum computer. Then entanglement between 
Alice and Bob does not make sense any more. 

An alternative (in the communication picture) is to let Alice and Bob have access to a noiseless 
channel, rather than entanglement. This channel serves as a catalyst and is returned at the end 
of the protocol. An [[n, k;c]] catalytic quantum error correcting (CQEC) code is defined by 

• An encoding isometry £ : C s>k <g> £®° — > C® n 

• A decoding CPTP map V : C® n ® C® c -> C® k 
such that 

V o {Af <g> id® c ) og = id® k = id® fc " c id® c . (39) 
The above may be written as a resource inequality 

(AO +c[q^q\>(k-c)[q^q\+c[q^q\. (40) 

The net rate of the CQEC code is given by (fc — c)/n. Figure [S] shows how any [[n, k; c]] EAQEC 
code {£,T>) gives rise to a CQEC code with the same parameters. This construction may be 
understood in terms of resource inequalities. The simple protocol called entanglement distribution 
written as 

c[q-> q] > c[qq], 

creates c ebits of entanglement by sending half of a locally prepared state |$)® c through the channel 
id® . The CQEC code is obtained by combining entanglement distribution with the EAQEC code: 



(AO + c [g -> q] 
>(JV)+c[qq] 
>k[q^q] 

= (k - c)[q -> q] + c[q -> q]. 

Assume now that Alice and Bob have access to m copies of the chanel Af. Performing the 
CQEC protocol m times in parallel (i.e. using the code (£® m , £>® m )) gives 

m(Af) + mc [q — > q] > m(k — c) [q — > q] + mc [q — > q] . 

The size of the catalyst can actually be reduced from mc to c: 

m(Af) + c [q —> q] > m(k - c) [q — ► q] + c [q —> q]. (42) 
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Figure 5: A catalytic quantum error-correcting code. 
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Figure 6: The inductive step. 



The proof is by induction. The statement is trivial for m = 1. For the inductive step, assume true 
for to. Then l|42|l holds for m + 1: 



(m+l)(jV) +c[q^q] 
= (A/) +m(A0 + c[q -> q] 

> (M) + m{k - c) [q -> g] + c [q -> q] 

> m(k — c) [q — > g] + (fe - c) [g — > g] + c [g 



(43) 



g]. 



A more conventional formulation of this catalyst reduction is given in the lemma below. 

Lemma 4.1 If \3y\) is satisfied then for any non-negative integer to there exists a CQEC code 
(£ m ,T> m ) for the channel M® m in the sense that 



T> ri 



(JV 6 



i id® c ) o £„ 



.id 6 



□ 



Proof The inductive step is shown in the Figure EJ 

The above construction is rather sensitive to perturbations. If in any particular block a channel 
worse than Af is experienced, the resulting channel will not be pure and the next block will start 
with an impure catalyst. 

One may rightly ask about where one could obtain a catalyst to begin with. A possibility is 
to use an ordinary c = QEC code. This is shown in Figure An [n, k; c] CQEC code for the 
channel M combined with a [n 1 , c; 0] QEC code for the channel N 1 gives an [n + n', k; 0] QEC code 
for the channel A/"® J\f'. The combined code can be used as a catalyst for an even larger code. In 
this way a sizeable catalyst can be built up pretty quickly. 
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Figure 7: Constructing a QEC code from a seed QEC code and a CQEC code. 



It is worth looking at this construction from a purely mathematical point of view. Let C C 
(Z2) 2 " and C C (Z2) 2 ™ be the symplectic codes corresponding to the [n, k; c) CQEC code and 
[n',c;0] QEC code, respectively. Let H and H' be the respective parity check matrices, as in 
section l2~!Tl Note that C' ± = iso(C"" L ). Let Ui, i = 1, . . . , c, be vectors in (Z2) 2 ™ which, together 
with a basis for C , form a maximal n'-dimensional isotropic subspace of (Z2) 2 ™ . Recall the 
notation = (ej|0) € (Z2) 2c . Let T be a symplectomorphism such that T(gj) = u^. Define the 
(n — k + c) x 2n' matrix B' = T(B) with B defined as in (|27|l and £ = n — k — c. Note that the 



rows of B' are in C . Then 



B' 

w 



H, 

c) x2n 1 

is the parity check matrix for the combined [n + n', k; 0] QEC code. By construction, is must be 
dual-containing. 



5 A variation on EAQEC codes 

One lesson learned from quantum Shannon theory |KS| is that catalytic and non-catalytic codes 
have similar performance. In this section we mimic the quantum Shannon theoretical construction 
from First we construct codes for sending classical information with entanglement assistance. 
Then we make these protocols coherent in the sense of ^3 EH to obtain a variation on EAQEC 
codes in which entanglement is generated as well as quantum communication. The end result is 
what we will call "type II" EAQEC codes, which can be constructed without the machinery of 
symplectic linear algebra. 

5.1 EA-codes for sending classical information 

The communication scenario again involves two spatially separated parties, Alice and Bob. The 
resources at their disposal are a noisy channel Af : C® n — ► C® 11 and the shared c ebit state l^)® 13 . 
Now Alice wishes to convey an element of (F4) fc perfectly to Bob using the above resources. A 
protocol which does this is called a [n, k; c]^ entanglement- assisted classical error correcting code, 
or EACEC code for short. We write the above as a resource inequality 

(Af) +c[qq]>2k[c->c]. (44) 

The factor of 2 accounts for the conversion from quaternary to binary. 

Recall the isomorphism 7 : (F4)™ — ► (Z2) 2 ". It allows us to, with a slight abuse of notation, 
speak of error sets S C (F4) ra , and Pauli matrices A^ a , a S (F4)". Let S C (F4)™ be the support of 
Af. An easy modification of Theorem 12.11 ensures that correctly decoding the message for the set 
of channels {N a : a 6 S} suffices for the correct decoding of Af. The notion of distance for EACEC 
codes is equivalent to the one for classical quaternary codes. An [n, k; c] 4 EACEC code of distance 
d is called a [n, k, d; c\± EACEC code. 

Proposition 5.1 If there exists an [n, k]i classical code (over F4J which corrects the error set 
S C (F4)™, then there exists an [n, fc;n]4 EACEC code which corrects the same error set. 
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Figure 8: Reduction from an EACEC code to a classical code over ¥4 



Proof We will show that superdense coding establishes an equivalence between a quantum Pauli 
error A c and a classical error c. 

Assume c = 0, corresponding to no error. Alice superdense encodes b by performing N\, 
on her half of |$)®". Bob performs a measurement in the { | < &t.)( < 3?b | : b £ (H^) 71 } basis, where 
|$ b ) = (A b ® thus decoding b. 

If the channel is N c for some c £ S, then Alice's effective encoding becomes A c Ab which is a 
representative of [Ab+ C ]- Bob's measurement will reveal b + c instead of b. This is the message 
with a classical error c £ S. The encoding preparation, followed by quantum error N c and decoding 
measurement, simulates the noisy classical channel b 1— * b + c. The theorem now follows, since the 
classical code can correct any error c £ S . □ 

Thus there is a direct correspondence between [n, k, d\i classical codes and [n,k,d;n] EACEC 
codes. On the other hand, in Section|2|we saw that an [n, k, d\± classical code defines an [[n, 2k — 
n + c, d; c]] EAQEC code. In the next subsection we show how to construct a variation on an 
[[n, 2k — n + c, d: c]] EAQEC code from an [n, k, d; n]^ EACEC code via "coherification." 



5.2 Coherent EACEC codes 

At this point we need to introduce one more resource, coherent communication |19| . Let {|0), |1)} 
denote a preferred basis for a qubit system. The isometric channel which implements the change 
of basis 

A 2 :\i) A ^\i) A \i) B , i = 0,l 

is called the coherent bit (or cobit) channel. The superscript A denotes a system held by Alice and 
B denotes a system held by Bob. It is regarded as a coherent version of a classical bit channel. 
Viewing it as a resource, we use the symbol [q — > qq\. Coherifying a protocol is a broad notion 
marked by replacing classical communication by coherent communication [I'M I19j . It was shown 
in |19j that superdense coding can be made coherent, i.e. that the following resource inequality 
holds: 

[q^q] + [qq] >2[q^qq\. (45) 

Consider an [n, k, d; n] EACEC code, given by l]44[l. It can also be made coherent thanks to its 
connection to superdense coding. In other words, l|44J) can be upgraded to 

(AO +n[qq] > 2k[q -> qq}. (46) 

An explicit circuit implementing this resource inequality is given in Figure H The states 
{|a) : a <E (F4) fe } form a basis for a 2k qubit space. {N c } is a Pauli matrix whose index c £ (F4)™ 
is in the support of Af. H4 is the (n ~ k) x n quaternary parity check matrix for the classical 
[n, k,d]i code which corrects all such c. G4 is the corresponding n x k generator matrix such that 
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coherent preparation 

|o> 






\H 4 c) 






|G 4 a + c)_ 


u 








\a) 





coherent measurement 



V 



\H A c) 



u 



\a) 



Figure 9: The circuit implementing a coherent EACEC code. The upper right hand corner defines 
U in terms of the quaternary code with parity check matrix H4 and generating matrix G4. 



H4G4 = 0( n _fc)xfc- The box in the upper right hand corner defines the 4™ x 4" unitary matrix 
U. There G4& S (F4) n is an encoded element a of (F4) fc . The unitaries V and W are given by 
V = E j evJMj\®N j and 

W(\cp)\0))= £<* + |(JVj<8>J)|¥>) 

Harrow |19| also exhibited a coherent version of quantum teleportation 1 , written as 

2k [q -> qq] + k [qq] > k [q q] + 2k [qq\. (47) 

Figure ITHl depicts a circuit implementing this resource inequality. 
Combining l|4t)|) with 14711 gives 

(TV) + (n + k) [qq] > k [q -> q] + 2k [qq]. (48) 

This differs from a hypothetical [[n, k;n — k]] EAQEC code given by 

(W) + (n — k) [qq] >k[q^q] 

in that an extra 2k [qq] is needed as a catalyst. We call this a type II [[n, k;n — k; 2k]] EAQEC 
code, and will refer to the EAQEC codes from Section 2 as type I EAQEC codes. A EAQEC 
code is not as versatile as regular EAQEC codes. The catalyst does not allow it to be converted 
into a catalytic QEC code, for example. Also, type II EAQEC codes appear to be limited to F4 
construction. 



This EAQEC code has the maximum value of c = n — k. 
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Figure 10: The circuit implementing coherent teleportation. 



As in the original Shannon theoretical result ^31 (Figure 0), type II EAQEC codes (|48(l can be 
combined with superdense coding 1)18)1 to give a catalytic version of an EACEC code 1)44)1 : 

(A/) + n[qq] + k[qq]>2k[c^> c]+k[qq]. 

This does not hold for type I EAQEC codes of Section[21 unless c equals its maximal value of n—k. 

6 Bounds on performance 

In this section we shall see that the performance of EAQEC codes is comparable to the performance 
of QEC codes (which are a special case of EAQEC codes). 

The two most important outer bounds for QEC codes are the quantum Singleton bound [211129) 
and the quantum Hamming bound 16 . Given an [[n, k,d\] QEC code (which is an [[n, k, d; 0]] 
EAQEC code), the quantum Singleton bound reads 

n — k > 2(d — 1). 

The quantum Hamming bound holds only for non-degenerate codes and reads 

y( n ) <2"- fc . 



3=0 



The proofs of these bounds [20 Ell are easily adapted to EAQEC codes with k — c playing the role 
of k. This was first noted by Bowen 5 in the case of the quantum Hamming bound, which now 
reads 

£ 3^ ( n \ < 2«- fe+c . 
j=0 ^ 3 ' 

The new quantum Singleton bound is 

n-k + c> 2{d - 1). 

Note that the F4 construction connects the quantum Singleton bound to the classical Singleton 
bound n — k > d — 1. An [n,k, d]i code saturating the classical Singleton bound implies an 
[[n, 2k — n + c, d; c]] EAQEC code saturating the quantum Singleton bound. 
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It is instructive to examine the asymptotic performance of quantum codes on a particular 
channel. A popular choice is the tensor power channel Af® n , where AT is the depolarizing channel 
with Kraus operators {-y/pol, ^/piX, ^Jp^Y, ^/p^Z}, for some probability vector p = (ft),Pi,P2,f?3)- 

It is well known that the maximal transmission rate R = k/n achievable by a non-degenerate 
QEC code (in the sense of vanishing error for large n on the channel AT® 71 ) is equal to the hashing 
bound R = 1 — H(p). Here H(p) is the Shannon entropy of the probability distribution p. This 
bound is attained by picking a random dual-containing code. However no explicit constructions 
are known which achieve this bound. 

Interestingly the F4 construction also connects the hashing bound to the Shannon bound for 
quaternary channels. Consider the quaternary channel a a + c, where c takes on values 0, u>, 1,57, 
with respective probabilities po,Pi,P2,P3- The maximal achievable rate R = k/n for this channel 
was proved by Shannon to equal R = 2 — H(p). An [n, k]i code saturating the Shannon bound 
implies an [[n, 2k — n + c; c]] CQEC code achieving the hashing bound! (Recall that the net rate 
of an [[n, k;c]\ CQEC code is defined as R — (k — c)/n.) Efficiently decodable modern classical 
codes, such as low-density parity check (LDPC) codes and turbo codes, are known to come very 
close to the Shannon bound. These codes are not at all guaranteed to be dual-containing |26|. 
which means that no hashing bound attaining QEC code can be constructed from them. However, 
a catalytic QEC code can. Directly investigating the symplectic structure of such modern codes 
will reveal the size of the catalyst. We expect that bootstrapping the method from Figure [7\ will 
enable us to construct a QEC code with similar properties. 



7 The [[3,1,3;2]] EAQEC code 

In this section, we will construct a [[3, 1, 3; 2]] EAQEC code and relate this code to Bowen's earlier 
result [5]- Consider the classical [3, 1, 3] quaternary code with parity check matrix 



Then 



Ha 



H = j(H 4 ) 





1 













(I 





1 


) 










1 














\ 


1 





1 






















1 


1 







V 








1 





1 





Following the proof of Theorem 1.1, we have 



ui = (1 1 0) 

u 2 = (0 1 1 0) 

U3 = (1 1 1 0) 

vi = (0 1 1) 

v 2 = (1 1 0) 

v 3 = (0 1 1 1), 



(49) 



(50) 



(51) 



and the hyperbolic pairs (ui,vi) and (ua,V2) span the rowspace of H. The simultaneous +1 
eigenstate of the commuting operators N Ui , i = 1,2, 3, is 

|000) = -±=(|000) + 1 110)). 
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Then 



000) 
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= N V1 1000) = 


010) 


= iV V2 1000) = 


Oil) 


= iv Vl+V2 |obo) 


100) 


= iV V3 |000) = 


101) 


= JV V1+V3 |000) 


no) 


-iv V2+V3 |oo6) 


m) 


= A f Vi+V 2 +V 3 K 



1 

71 
i 



V2 



(|101) -|011)) 



1 



v/2 



(-|101) + |011)) 



V2 



(|1H) + |001» 



1 

71 
i 

= 71 
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The encoding unitary Ur is therefore 



The logical and 1 codewords are 



( 1 





1 































1 





1 






















1 





1 







1 





1 































1 





-1 







1 





-1 
















1 





-1 






























1 





-1 





) 



\0 L ) = C/ T |0)|$+)® 2 = ~ ( |000)|00) + |001)|01) + |010)|10) + |011)|11) 



|1 L ) = U T \m+) m = \ ( |100)|00) + |101)|01) + |110)|10) + |111)|11) 



(52) 



(53) 



(54) 



Bowen's code [5] can be obtained by applying the following unitary to the codewords given above 

( 1 -1 -1 -1 -1 1 1 1 \ 

-1-1-11 1 1 1-1 
1-11 1-11-1-1 
-1-11-11 1-11 

1 -1 -1 -1 1 -1 -1 -1 

1 1 1-11 1 1-1 

-1 1 -1 -1 -1 1 -1 -1 

y —1 — 1 1 -1 -1 -1 1 -I J 



U B = 



2a/2 



(55) 



8 Table of codes 

In [Hj a table of best known QEC codes was given. Below we show an updated table which includes 
EAQEC codes. 
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The entries with an asterisk mark the improvements over the table from All these are 
obtained from Proposition 3.1. The corresponding classical quaternary code is available online at 
http : / / www . win . tue . nl/^aeb/voorlincod . html 

The general methods from [H] for constructing new codes from old also apply here. Moreover, 
new constructions are possible since the dual-containing condition is lifted. An example is given 
by the following Theorem. 

Theorem 8.1 a) Suppose an [[n,k,d;c]] code exists, then an [[n + 1, k — c + c' — 1,g?';c']] code 
exists for some c' and d' > d; b) Suppose a non- degenerate [[n, k, d;c]] code exists, then an [[n — 
1, k — c + c' + 1, d — 1; c']] code exists for some c' . 

Proof a) Let H be the {n — k + c x 2n) parity check matrix of the [[n, k, d; c]] code. The parity 
check matrix of the new [[n + 1, k — c + c' — 1, d'; c']\ is then 



H' 



( • • • 

1 • • • 1 1 
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1... 1 1 \ 

••• 




(56) 
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This corresponds to the classical construction of adding a parity check at the end of the codeword 
|27| . The additional rows ensure that errors involving the last qubit are detected. Sometimes the 
distance actually increases: for instance, the [[8,0,4;0]] code is obtained from the [[7, 1,3;0]] code 
in this way. 

b) We mimic the classical "puncturing" method £ 27. . Let C be the (n + k — c)-dimensional 
subspace of (Z 2 ) 2n corresponding to the [[n, k,d;c]} EAQEC code. Puncturing C by deleting the 
first Z and X coordinate, we obtain a new "code" C which is an (n + k — c)-dimensional subspace 
of (Z2) 2 (™ -1 - ) . This corresponds to an [[n — 1, A; — c+c' + 1, d— 1; c']] EAQEC code, as the minimum 
distance between the "codewords" of C decreases by at most 1. □ 



9 Discussion 

Motivated by recent developments in quantum Shannon theory, we have introduced a generalization 
of the stabilizer formalism to the setting in which the encoder Alice and decoder Bob pre-share 
entanglement (EAQEC codes). We have traced the male side of family tree of quantum Shannon 
theory, from EAQEC codes (corresponding to the father protocol) to catalytic quantum codes 
(corresponding to the quantum capacity) and EACEC codes (corresponding to the classical EA- 
capacity). Moreover, EACEC codes can be made coherent, providing an alternative to the EAQEC 
construction from Section 2. The most obvious question is whether we can do the same for 
the female side of the family tree ^21- Preliminary results [2H give a positive answer to this 
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question: entanglement distillation protocols assisted by quantum and classical communication 
can be constructed based on non-orthogonal symplectic codes. 

There are two practical advantages of EAQEC codes over standard QEC codes: 

1. They are much easier to construct from classical codes because they are not required to be 
dual-containing. This allows us to import the classical theory of error correction wholesale, 
including capacity-achieving modern codes. We plan to examine the performance of classical 
LDPC codes and turbo codes in the context of the catalyst size for EAQEC codes. 

2. The entanglement used in the protocol is a strictly weaker resource than quantum communi- 
cation. Thus comparing [[n, k, d; c]] EAQEC codes to [[n, k — c,d; 0]] QEC codes is not being 
entirely fair to former. The pre-shared entanglement could have been obtained from a two 
way entanglement distillation protocol which achieve higher rates compared than one-way 
schemes. In this sense, a large value of the catalyst c is viewed as advantageous, as it implies 
a higher qubit channel yield. 

If one is interested in applications to fault tolerant quantum computation, where the resource 
of entanglement is meaningless, high values of c are unwelcome because they require a long seed 
QEC code. We expect this obstacle to be overcome by bootstrapping. 

Another fruitful line of investigation connects to quantum cryptography. Quantum crypto- 
graphic protocols, such as BB84, are intimately related to CSS QEC codes. In [23 it is shown 
that EAQEC analogues of CSS codes give rise to key expansion protocols which do not rely on the 
existence of long dual-containing codes. 

We thank Graeme Smith for pointing us to references and ^5] • TAB acknowledges financial 
support from NSF Grant No. CCF-0448658, and TAB and MHH both received support from NSF 
Grant No. ECS-0507270. ID and MHH acknowledge financial support from NSF Grant No. CCF- 
0524811 and NSF Grant No. CCF-0545845. 
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